Pure Mott phases in confined ultra-cold atomic systems 
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We propose a novel scheme for confining atoms to optical lattices by engineering a spatially- 
inhomogeneous hopping matrix element in the Hubbard-model (HM) description, a situation we 
term off-diagonal confinement (ODC). We show, via an exact numerical solution of the boson HM 
with ODC, that this scheme possesses distinct advantages over the conventional method of confining 
atoms using an additional trapping potential, including incompressible Mott phases at commensu- 
rate filling and a phase diagram that is similar to the uniform HM. The experimental implementation 
of ODC will thus allow a more faithful realization of correlated phases in cold atom experiments. 
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In recent years there has been a burgeoning interest 
in the boson Hubbard (BH) model [TJE], spurred by the 
expectation [3] that modern laser cooling and trapping 
techniques lead to the realization of this model in cold 
atom experiments. Several experimental groups [4HTT] 
have achieved the BH model with bosonic atoms in op- 
tical lattices, observing evidence of the predicted Mott 
insulating and superfluid phases. Experimental probes 
include the condensate fraction and momentum distribu- 
tion [8], noise correlations [9] and RF spectroscopy [7]. 

The ground state phase diagram for the uniform BH 
model with chemical potential /i, onsite repulsion J7, and 
hopping parameter t is well known (Fig. [T]), and con- 
sists of lobes of incompressible Mott insulating phase (at 
commensurate fillings) surrounded by a regime of super- 
fluid order, with these states separated by quantum phase 
transitions [TJ [2] . However, present-day cold atom exper- 
iments always occur in the nonuniform environment of a 
smoothly- varying harmonic trap that confines the atoms 
to the optical lattice, so that the actual Hamiltonian con- 
sists of the BH model plus a parabolic single-particle po- 
tential. Early theoretical work [12] showed that this kind 
of diagonal confinement (DC) does not entirely reproduce 
the physics of the BH model. In particular, the model 
with DC is characterized by the absence of energy gaps 
at commensurate fillings, as opposed to the BH model. 
This results in the absence of true incompressible Mott 
insulating phases and of a real superfluid to Mott tran- 
sition in present-day experiments. 

The physical origin of the lack of incompressible Mott 
phases in the BH model with DC [3] [12] can be under- 
stood within the local density approximation (LDA). Ac- 
cording to the LDA, the harmonic DC potential is treated 
by replacing \i — >> \i — ^mQ 2 r 2 with r the distance from 
the center of the lattice, Q the trapping frequency and m 
the atom mass. With increasing radius, then, the system 
traverses the uniform-case phase diagram at fixed t and 
a spatially- varying \i (vertical arrow on Fig.0. Then, 
even if a system is locally in a Mott phase at the cen- 
ter of the trap, it must enter a superfluid phase with 



increasing radius (yielding the well-known shell struc- 
ture |6j [3 [I2j [13] ) , thereby being globally compressible. 
Recent work has used the LDA to infer the homogeneous 




FIG. 1: (Color online) Schematic ground state phase diagram 
of the uniform boson Hubbard (BH) model [T]. The arrows 
show two possible sequences of local phases encountered in an 
atomic system as the radius r from the center of the lattice 
increases. Diagonal confinement (DC) is created by a decreas- 
ing chemical potential /x, leading to a sequence of Mott and 
superfluid phases. Off-diagonal confinement (ODC) is created 
by a decreasing hopping parameter t, and can maintain the 
whole system in the same Mott lobe. Inset: A decreasing hop- 
ping parameter t can be obtained by increasing the intensity 
of the optical lattice as one moves away from the center. 

BH phase diagram from the DC experimental results [TT] 
and also to place a given DC system on the phase dia- 
gram of the uniform BH case [14]. Nevertheless, it is still 
of great interest to establish truly incompressible Mott 
phases of bosons in cold atomic gas experiments. 

In this Letter, we propose that truly incompressible 
Mott phases can be achieved via off-diagonal confinement 
(ODC) of the bosons where the cloud is confined via a 
spatially varying tunneling matrix element t(r), which 
vanishes at the boundary of the system. This can be re- 
alized, for example, by having the optical lattice depth 
increase with increasing radius, as illustrated in the inset 
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of Fig. [T] Although a system with ODC is still spatially 
inhomogeneous, it can exhibit a true gapped and incom- 
pressible Mott phase at commensurate fillings. This can 
be easily understood by extending the above LDA argu- 
ment to a spatially- varying t(r) instead of ji{r). Thus, 
within the LDA, a vanishing t(r) with increasing r corre- 
sponds to traversing the uniform-case phase diagram at 
fixed chemical potential. In particular, a system that is 
locally in a Mott phase at the center of the lattice will 
remain in the same Mott phase for all radii up to the 
boundary where t(r) becomes zero. 

Thus the ODC setup reproduces important aspects of 
the BH model, namely the presence of gaps and true 
Mott phases at commensurate fillings, as we show via an 
exact numerical calculation for a one-dimensional model. 
In addition, superfluid phases with ODC are found to 
lead to condensates (quasicondensates in one dimension) 
with greater population than the model with DC, and 
therefore can realize clean condensates. 

Let us briefly discuss the experimental viability 
of achieving off-diagonal confinement via a spatially- 
dependent t. A uniform periodic lattice of cold atoms 
yields a tunneling matrix element [15] 
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with E r the recoil energy (E r 

lattice wave vector, m the atom mass) and Vo the opti- 
cal lattice depth. The exponential sensitivity of t on Vo 
implies that even a modest increase in Vo with increas- 
ing radius will yield a rapid suppression of t, confining 
the cloud to the optical lattice. Taking parameters from 
Ref. [H which has E r ~ 2.13 x 10" 13 J, this sensitivity 
is illustrated by noting that the hopping time t ~ h/t 
changes by two orders of magnitude as Vo/E r is varied 
(from 0.3ms at V /E r = 3 to 26ms at V /E r = 22). 

Although conventional optical lattices created with in- 
terfering laser beams typically possess a uniform t, re- 
cently the Greiner group has developed a novel holo- 
graphic method that can produce an arbitrary optical po- 
tential using a mask [16] . An optical lattice that produces 
a spatially-varying t will also yield a spatially-varying 
Coulomb repulsion U. However, examining the expres- 
sion for this quantity [15] U = ^8/7rkaE r (Vo/E r ) 3 / 4 , 
with a the atom scattering length, we note that U is only 
power-law sensitive to Vo. Therefore, a spatially- varying 
Vo will lead to relatively small variations of U and larger 
changes in t. Moreover, variations of U can be taken 
into account by generalizing the LDA picture to a vary- 
ing t/U, yielding same conclusions. Thus, we expect that 
such off-diagonal confinement of atoms in an optical lat- 
tice will be experimentally realizable in the near future, 
adding to the rapidly-growing cold-atom toolbox [TTj . 

The above LDA argument suggests that experiments 
using ODC will achieve a true Mott phase for any 



spatially- varying t(r). To illustrate this, we consider the 
following one-dimensional model: 



H 



Uj (a\aj + a)ai) + - ^ hi (hi - l) (2) 




FIG. 2: (Color online) The density as a function of the chem- 
ical potential, for [7 = 8. The ODC model exhibits Mott gaps 
at commensurate fillings as in the periodic case (inset). Re- 
sults for the DC model with a trapping parameter W = 0.008 
are shown for comparison. 

The sum over i runs in the range [0; L — l] , with 
L the number of lattice sites, and the sum over Uj} 
is restricted to pairs of first neighbors. The bosonic 
operators a\ and a i create and annihilate a parti- 
cle on site z, respectively, and hi = a\a i . The lo- 
cal hopping parameter tij follows an inverted parabola, 
tij = t(i + j + 1)(2L — i — j — 1)/L 2 , with a maximum 
value of t at the center of the lattice, tz,/2-i,L/2 = t, 
and a vanishing value at the edges (£-i,o = = 0). 

Finally, we choose t = 1 in order to set the energy scale. 

We solve the model exactly by means of quantum 
Monte Carlo simulations, by making use of the Stochas- 
tic Green Function (SGF) algorithm [18 with directed 
updates [19]. The SGF algorithm can be implemented in 
the canonical ensemble as well as in the grand-canonical 
ensemble. In the following, we take advantage of this fea- 
ture. When working in the canonical ensemble, the num- 
ber of particles N is chosen and remains strictly constant 
during the entire simulation. In the grand-canonical en- 
semble, we add the usual term —fiN to the Hamilto- 
nian, with N = J^fti, in order to control the mean 
number of particles N = (iV) via the chemical poten- 
tial ji. In the following we consider a lattice with L = 70 
sites, since this is the typical size currently accessible 
in one-dimensional optical lattice experiments. We have 
used increasing values of the inverse temperature /?, and 
found that the ground-state properties are captured with 
f) = 20. All presented results correspond to this inverse 
temperature. We begin by investigating the behavior of 
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the density p = N/L as a function of the chemical po- 
tential p. Here it is convenient to make use of the grand- 
canonical ensemble, since the chemical potential is the 
control parameter. Fig. [2] displays results for U = 8. 
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FIG. 3: (Color online) The local density and local compress- 
ibility, for U = 8. A superfluid cloud grows in the center of 
the lattice, as the number of particles increases (Panel A). 
Contrary to the DC model, Mott regions in the ODC model 
form at the edges of the lattice (Panel B), then extend towards 
the center. At commensurate filling (Panel C), the phase is a 
pure incompressible Mott insulator. Adding particles above 
commensurate filling breaks the Mott plateau in the center, 
which becomes locally superfluid (Panel D). 

A comparison between the results in the main panel 
and the inset in Fig. [2] shows that while the usual model 
with DC does not display any incompressible phase, the 
model with ODC has pure Mott gaps at commensurate 
fillings as in the homogeneous case. For these fillings, 
the number of particles per site sticks to integer values 
and is constant over the whole lattice. Note also that 
the gaps of the ODC model and homogeneous BH model 
are similar in size. The results for the DC model corre- 
spond to a parabolic potential — L/2) 2 fii with 
W = 0.008. For this case, the size of the lattice has been 
increased to L = 100 in order to avoid an overflow of the 
particles. Note that such an overflow never occurs in the 
ODC model, since the local hopping parameter vanishes 
at the edges of the lattice. 

Fig. [3] shows typical density profiles and the associated 
local compressibility, defined as the response of the local 
number of particles to a change of the chemical potential 

CSII, «i = ^ = f [(ni(r)N) - (ni(T))(N)]dT. For 
oi 1 Jo 

pure Mott phases this quantity vanishes at commensurate 
fillings for each lattice site i. 

The density profiles in Fig. [3] can be qualitatively un- 
derstood within the above LDA argument, in which in- 
creasing radius corresponds to a decreasing t and a hor- 
izontal path in the homogeneous BH phase diagram of 
Fig. [I] Thus, Fig. [3]A shows a system that is locally 
superfluid in the center, and enters the vacuum as the 



radius increases and the ratio t/U decreases. For greater 
/i, Fig. corresponds to a path beginning in the super- 
fluid phase and crossing into the n = 1 Mott phase of 
Fig. [l] therefore the system is still globally compressible 
with a local Mott phase at the edge. The incompress- 
ible Mott phase in Fig. |3p corresponds to a LDA path 
that remains entirely inside one Mott lobe. Fig. [3p at 
higher p is qualitatively similar to Fig. [3J3, but with the 
inner superfluid having a higher local density than the 
Mott region at the edges. Finally, all density profiles 
corresponding to the p values explored in Fig. [2] are rep- 
resented on Fig. [4] The figure clearly shows the Mott 
plateaux at p = 1 and p = 2 extending over the whole 
lattice. 




FIG. 4: (Color online) The density as a function of the posi- 
tion and the chemical potential for U — 8. The Mott plateaux 
at p — 1 and p — 2 appear clearly, and extend over the whole 
lattice. 
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FIG. 5: (Color online) The ground state phase diagram of the 
ODC model given by Eq. pi, showing lobes of incompressible 
pure Mott phase (like Fig~[3p), and regions of compressible 
superfluid and compressible coexisting Mott and superfluid 
phases (like Figs. [3}\, [3^, and[3p). 

Next, we work in the canonical ensemble at fixed den- 
sities and scan the values of the onsite repulsion param- 
eter U for which the system exhibits the Mott insulating 
state. In the canonical ensemble the ground state chemi- 
cal potential is not a control parameter, but instead it is 
computed as p(N) = E(N + 1)-E(N). Fig. [5] shows the 
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phase diagram which exhibits pure Mott lobes similar to 
those of the homogeneous BH model [2] , which again are 
not present in the DC model [12]. In addition, new lobes 
where the phase is purely superfluid arise as the onsite 
interaction decreases. Another particularity of the ODC 
model is that coexistence of Mott and superfluid can oc- 
cur at any value of the onsite repulsion U . This mixed 
state is absent in the uniform BH model; and in the DC 
model there is a critical value of U [21] below which the 
system is purely superfluid. 

Finally, we make connection with an experimental 
probe by plotting in Fig. [6] the maximum intensity I m ax 
of the interference pattern (the population of the quasi- 
condensate) as a function of the chemical potential. The 

Visibility V = {I max ~ ^min) / '(Imax + Imin) is also shown. 

The variations in the visibility are stronger with ODC 
than with DC, since in the former case we are in the 
presence of a true superfluid-to-Mott transition. The 
maximum intensity curve also reveals another interest- 
ing feature of the ODC model: condensates have greater 
population than the DC model. Thus, the experimental 
realization of the ODC model should improve the quality 
of condensates. 




FIG. 6: (Color online) The maximum intensity of the inter- 
ference pattern and the visibility as functions of the chemical 
potential, for [7 = 8. 



Mott phase will pave the way to further experimental 
simulations of the Mott-superfluid transition, a problem 
of broad interest due to its connections to condensed- 
matter systems such as the high-temperature supercon- 
ductors [22] . In addition, our simulations show that the 
model leads to highly populated condensates, which are 
of great interest in the production of coherent matter 
waves. 

To achieve the ODC, it may be easiest to start with a 
conventional DC cloud (a system that has already been 
taken to low temperatures) and adiabatically apply the 
ODC optical potential. We have studied the entropy of 
the DC and ODC phases in a hardcore-boson model at 
finite temperatures, finding, for comparable system size, 
a comparable entropy in the ODC and DC states, im- 
plying that a low temperature can be maintained upon 
adiabatically moving from DC to ODC. Additionally, one 
may ask how equilibrium is established with ODC given 
that the local hopping time r ~ %/t is spatially varying; 
rapid equilibration may be stymied if much of the system 
possesses a large value of r. To address this it may be 
necessary to engineer ODC to have Vq increase rapidly 
over a few lattice sites near the edge. Given the above- 
mentioned rapid variation of r with Vq, this should be 
feasible. 

The flexibility of the techniques used in present-day ex- 
periments allows the realization of our model, and should 
open new perspectives on the physics of Mott phases. 
Future work will utilize the Stochastic Green Function 
(SGF) algorithm [18] [19] to compute other experimen- 
tally relevant quantities such as the momentum distribu- 
tion [8 , noise correlations [9] and RF spectroscopies [7]. 
Finally, we note that our model can also be easily gener- 
alized to fermionic systems where we expect a similarly 
rich behavior. 
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In this Letter, we proposed a novel scheme for confin- 
ing atoms on optical lattices based on off-diagonal con- 
finement, in which the Bose-Hubbard tunneling matrix 
element is engineered to vanish at the system bound- 
ary. We show by using an exact quantum Monte Carlo 
method that our model reproduces important aspects of 
the boson Hubbard model. In particular, truly incom- 
pressible bosonic Mott phases occur at commensurate 
fillings and strong interactions, an important require- 
ment for quantum-computing [17 . Indeed, optical lat- 
tices provide intrinsically scalable and well-defined sets 
of qubits. These qubits must be initialized in a "zero" 
state, which corresponds to having one particle per site. 
This is precisely achieved by having the system in a pure 
Mott phase. Additionally the establishment of a global 
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